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ABSTRACT 

Variable stars have been identified for the first time in Fornax 4, the globular 
cluster located near the center of the Fornax dwarf spheroidal galaxy. By applying 
the image subtraction technique to B,V time series photometry obtained with 
the MagIC camera of the 6.5-m Magellan/Clay telescope and with the wide field 
imager of the 4-m Blanco/CTIO telescope, we detected 27 RR Lyrae stars (22 
fundamental mode, 3 first overtone, and 2 double-mode pulsators) in a 2.4' x 2.4' 
area centered on Fornax 4. The average and minimum periods of the ab-type RR 
Lyrae stars, (P a b) = 0.594 d and P a b,min = 0.5191 d, respectively, as well as the 
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revised position of the cluster in the horizontal branch type-metallicity plane, 
all consistently point to an Oosterhoff-intermediate status for the cluster, unlike 
what is seen for the vast majority of Galactic globular clusters, but in agreement 
with previous indications for the other globular clusters in Fornax. 

The average apparent magnitude of the RR Lyrae stars located within 30" 
from the cluster center is (V(RR)) = 21.43 ± 0.03 mag (a = 0.10 mag, average 
on 12 stars), leading to a true distance modulus of /i = 20.64 ± 0.09 mag or 
/i = 20.53 ± 0.09 mag, depending on whether a low ([Fe/H] = —2.0) or a 
moderately high ([Fe/H] = —1.5) metallicity is adopted. 

Subject headings: galaxies: dwarf — galaxies: individual (Fornax) — globular 
clusters: individual (Fornax 4) — stars: horizontal branch — stars: variables: 
other — techniques: photometry 

1. Introduction 

This is the first in a series of papers devoted to the detailed and comprehensive study of 
the variable star population in the field and globular clusters of the Fornax dwarf spheroidal 
galaxy (dSph) (The Fornax Project, Clementini et al. 2006), encompassing the galaxy's 
classical instability strip from the Anomalous Cepheid (V ~ 19 mag) to the Dwarf Cepheid 
(V ~ 25 mag) region (Poretti et al. 2006, 2007). Here, we report on results from the first 
variability study in Fornax 4 (For 4), the most compact and centrally located of the Fornax 
globular clusters (GC's), for which it has been suggested that it might actually be the nucleus 
of the Fornax dSph galaxy (Hardy 2002; Strader et al. 2003). In the following papers of this 
series, we will report on our study of the variable stars identified in the Fornax clusters #2, 
3 and 5, as well as in the survey of about a 1 deg 2 of the galaxy field. 

Fornax and Sagittarius are the only dSph's known to host globular clusters (GC's). The 
Fornax dSph GC system contains 5 GC's (Hodge 1961, 1965, 1969). Buonanno et al. (1998, 
1999) published color-magnitude diagrams (CMDs) of the Fornax GC's based on Hubble 
Space Telescope (HST) observations from which they estimated the clusters' relative ages and 
metallicities. They concluded that For 4 differs both in horizontal branch (HB) morphology 
and age from the other GC's in Fornax: it is a metal-poor cluster ([Fe/H] = —2.01 ±0.20, on 
the Zinn & West 1984 metallicity scale; Buonanno et al. 1999) with a red HB and seems to be 
3-4 Gyr younger than the other clusters. They also noted that the cluster CMD is virtually 
identical to that of the "outer halo" Galactic GC Ruprecht 106 (Rup 106), an unusual cluster 
often considered to have once been a member of a dwarf galaxy tidally disrupted by the Milky 
Way (e.g. Brown, Wallerstein, & Zucker 1997; Buonanno et al. 1999; Bellazzini, Ferraro, & 
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Ibata 2003). 

More recently, Strader et al. (2003) used low-resolution, integrated Keck spectra to 
measure metallicities and infer ages for the five GC's in Fornax. The new ages and metallic- 
ities are different than found by Buonanno et al. (1998, 1999). According to Strader et al., 
the metallicity of For 4 is [Fe/H] = -1.5 ± 0.12 (on the Zinn & West 1984 metallicity scale), 
and they estimate that the age of For 4 is the same as that of Fornax 1,2, and 3, all older 
than Fornax 5 by 2-3 Gyr. 

Because GC's are the oldest individual entities observable in galaxies, they provide 
invaluable insight into the formation history of the systems to which they belong. A puzzling 
feature of the Galactic GC's is that they sharply divide into two distinct types according 
to the mean periods of their RR Lyrae stars and the relative proportions of fundamental- 
mode (RRab) and first-overtone pulsators (RRc; Oosterhoff 1939). In the Milky Way (MW) 
Oosterhoff type I (Ool) clusters have (P a b) — 0.55 d, whereas Oosterhoff type II (OoII) 
ones have (P a b) — 0.66 d (Clement et al. 2001). Galactic GC's belonging to different 
Oosterhoff types may have different kinematical and spatial distributions, possibly resulting 
from different accretion/formation events in the halo (van den Bergh 1993; Yoon & Lee 
2002, Miceli et al. 2007). The existence of the Oo groups within the MW GC's can thus 
provide clues about the halo formation process. However, we do not know if the Oosterhoff 
dichotomy is a general characteristic of old stellar populations in galaxies, or if it is a peculiar 
phenomenon due to the particular evolutionary history of the MW. Studying the Oosterhoff 
properties of GC's belonging to the dwarf companions of the MW thus plays a crucial role 
in identifying the "building blocks" of the Galactic halo and in understanding which galaxy 
formation scenario, whether the fast monolithic free-fall collapse of protoclouds (Eggen, 
Lynden-Bell, & Sandage 1962; but see also Sandage 1990) or the hierarchical merging and 
continual accretion of lower-mass protogalactic fragments (Searle & Zinn 1978), is dominant 
(see, e.g., Catelan 2007). 

Mackey & Gilmore (2003a) identified candidate RR Lyrae stars in four of the Fornax 
dSph GC's (namely clusters # 1, 2, 3, and 5) based on archival HST observations. They de- 
termined periods by fitting template RR Lyrae light curves to their data and concluded that 
the Fornax clusters are unusual in that their RR Lyrae populations have mean characteristics 
intermediate between the Oosterhoff groups. 

For 4 has never been surveyed for variability, in spite of clear indications, from its HB 
morphology (Buonanno et al. 1999), that it should contain RR Lyrae stars. In this paper we 
present results from an extensive study of the variable star population in this extragalactic 
GC that resulted in the discovery of 29 variable stars, of which 27 are certainly RR Lyrae 
stars. The average luminosity of the RR Lyrae stars is used to measure the cluster distance. 
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The periods, amplitudes, and period-amplitude distributions allow us to define the OosterhofT 
type (Oo-type) of For 4 and to verify whether it conforms to the Oosterhoff dichotomy shown 
by the Galactic GC's. 

Observations and the adopted data reduction techniques are presented in Section 2. The 
variable star identification procedure is described in Section 3. The pulsation properties of 
the For 4 variable stars and their membership in the cluster are discussed in Section 4. In 
Section 5 we discuss the Oosterhoff classification and compare the For 4 Oo-type to that of 
the other GC's in Fornax and to the MW GC's. In Section 6 we derive the cluster distance 
from the average luminosity of its RR Lyrae stars. Final results are summarized in Section 7. 

2. Observations and data reductions 

B, V time series photometry of For 4 and of its surrounding field was obtained with the 
MagIC camera of the Magellan/Clay 6.5-m telescope on two nights in November 2003 and 
on two nights in December 2004. 

We covered a total field of view of 2.4' x 2.4', centered on the cluster at RA=2:40:07.33, 
DEC=— 34:32:17.4 (J2000). Nights were photometric with seeing conditions varying from 
0.45" to 0.65". We acquired a total number of 58 V and 19 B images. Average exposure times 
were 500 seconds for the V frames and 700s for the B ones. Observations of the standard 
fields T Phe and Ru 149 (Landolt 1992) were obtained on the same nights, to calibrate the 
data to the standard Johnson-Cousins photometric system. 

The Magellan data were complemented by observations of the Fornax dSph field con- 
taining For 4 obtained with the wide field imager of the 4-m Blanco/CTIO telescope. Data 
were collected over three runs for a total of 5 nights spread over three years, from 2003 to 
2005. The images were taken with the Mosaic II imager, which has a 36' x 36' field of view 
centered at RA=2:41:07.29 and DEC=-34:16:35.6 (J2000) with a plate-scale of 0.27"/pixel. 
The cluster was imaged on amplifier 1 of the 16-amplifier mosaic. The conditions during 
the three runs were photometric or nearly photometric. The seeing varied over the range of 
0.8" to 1.9". The exposure times varied from 100s to 600s in V and 200s to 600s in B with 
100s and 200s being typical, respectively. The CTIO dataset consists of 142 V and 61 B 
images. Logs of the observations and details of the instrumental set-up at the two telescopes 
are given in Table 1. 

The Magellan frames were bias-subtracted and flatfield-corrected using the MagIC tool 
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in irafQ They were then reduced with DAOPHOT-ALLSTAR. Time-series photometry 
for all the stars was produced with ALLFRAME (Stetson 1994, 1996). The photometric 
precision at the HB level is 0.01 mag. The photometric calibration was derived using the 
standard stars observed at the Magellan telescope on the night of Nov. 13, 2003 (UT). 

A set of linear calibration relations was computed: 



B = b' + 0.051(5-10 + 26.87, 
V = u'- 0.049(5- V) + 26.81, 

where b' and v' are the instrumental magnitudes normalized to 1-s exposure and corrected 
for atmospheric extinction by adopting the mean extinction coefficients ks = 0.22 and 
ky = 0.12. The zero point uncertainties of the calibration relations are estimated to be of 
the order of 0.05 mag in B and V, and 0.03 mag in (B — V), from the night-to-night scatter 
of the standard stars. These relations were used to calibrate the instrumental magnitudes, 
and aperture corrections were applied based on photometry of relatively isolated stars. 

The Magellan observations of For 4 produced a CMD reaching V ~ 25 mag. This CMD 
is shown in Figure 1. Although the cluster area observed with the Magellan telescope mainly 
contains stars within the For 4 tidal radius, the cluster CMD is clearly contaminated by stars 
belonging to the Fornax field, since the cluster lies in the center of the Fornax dSph galaxy. 
On the other hand, no reliable photometry could be obtained for stars inside a radius of 7" 
from the cluster center, due to crowding. The CMD in Figure 1 was drawn cutting off this 
central region. 

The images taken at the 4-m Blanco telescope were reduced using the SMSN pipeline 
(Rest et al. 2005) originally created for the SuperMACHO and ESSENCE projects. The 
SMSN pipeline reduces the raw images using the standard IRAF mosaic reduction package, 
MSCRED, as well as its own routines. As part of this reduction, the images were cross-talk 
corrected, astrometrically calibrated, and split into the 16 separate amplifiers. The cross-talk 
correction is necessary to remove ghost images of bright stars that can appear as a result 
of the electronics of the mosaic imager. Since the images are optically distorted due to the 
large field of view, they were then de-projected into a tangential plane using a customized 



1 IRAF is distributed by the National Optical Astronomical Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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version of Swarp (Bertin et al. 2002). DoPHOT (Schechter, Mateo, & Saha 2003) was then 
used on the frames to get profile- fitted photometry of the stars (and other objects), with 
photometric precision of 0.03-0.05 mag for stars on the HB. Local standards were used to 
transform individual images to the standard system. Final photometric calibration was then 
obtained by adjusting to the zero points of the Magellan dataset. 

3. Variable star identification 

Variable stars in For 4 were identified by applying the Image Subtraction Technique, 
as performed by ISIS2.1 (Alard 2000), to the 2003+2004 Magellan V dataset (58 frames 
in total), producing a catalog of about 400 candidate variable stars. The ISIS catalog was 
cross-identified against the DAOPHOT/ALLFRAME photometric catalog, leading to light 
curves in magnitude scale and confirming the variability for 20 stars. A further 9 variable 
stars were recovered in the cluster's central region where we lack reliable PSF photometry. 
For these stars we analyzed only the differential flux light curves produced by ISIS. 

Variable stars were independently identified from the CTIO dataset, using difference 
image analysis (Phillips & Davis 1995; Alard & Lupton 1998; Alard 2000). We used the 
HOTPANTS package (Becker et al. 2004) to match the PSF of the image and template by 
applying a spatially varying kernel to the image with the better seeing before the template is 
subtracted. One of the main problems with the image differencing approach is that there are 
more residuals, e.g. cosmic rays and bleeds, than genuinely variable objects in the difference 
image. Therefore standard profile-fitting software like DoPHOT has problems determining 
the proper psf used to perform photometry in the difference image. The difference image 
was analyzed with our customized version of DoPHOT, where we forced the PSF to be the 
one determined for the original, flattened image. Applying this a priori knowledge of the PSF 
helped to guard against bright false positives, such as cosmic rays and noise peaks, which 
generally do not have a stellar PSF. 

By matching the coordinates we were able to cross- identify 20 of the variable stars 
identified from the Magellan dataset. Fifteen of them had been independently recognized as 
variable stars from the CTIO data. Due to the high crowding no reliable match could be 
obtained for variable stars in the cluster central region. The time series B, V photometry of 
the 20 confirmed variable stars with light curves in magnitude scale is provided in Table 2 
and is available in electronic form in the electronic edition of the journal. 
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4. Period Search and Pulsation Properties of the Variable Stars 

The time series of the 400 candidate variable stars were analyzed independently by two of 
us (GC and EP) using two different codes and then compared. GC used GrATiS (Graphical 
Analyzer of Time Series), a custom package developed at the Bologna Observatory by P. 
Montegriffo (see Di Fabrizio 1999 and Clementini et al. 2000), which uses both the Lomb 
periodogram (Lomb 1976, Scargle 1982) and the best fit of the data with a truncated Fourier 
series (Barning 1963). 

EP used the iterative sine-wave least-squares method (Vanicek 1971), combined with a 
Fourier series decomposition limited to significant terms only (Poretti 2001). Visual inspec- 
tion of the intranight light curves was also performed. After having compared the results, 
only candidate variable stars confirmed by both methods have been accepted (29 stars in 
total). Our time series extends over three consecutive nights in 2003, 4 nights in 2004, and 
1 further night in 2005 (see Table 1). Variable stars with periods of several days or longer 
would not be detectable in the data set of each individual year, since the long period will 
show effects only on longer timescales. Such long-period variable stars, if present, would 
manifest themselves as stars having different magnitudes in the three years of our survey; 
since we have only 4-5 epochs, they would appear as periodic variable stars with a promi- 
nent low-frequency peak in the power spectra. Indeed, we do not have such candidates in 
our sample. At the same extent, non-periodic variable stars with timescale of 1-5 days will 
result in periodic variable stars, since any periodogram will find a (false) period when a 
non-periodic variable is sampled at 4-5 epochs only. Irregular, rapid variable stars are of 
course undetectable, but they are also not very probable. In conclusion, it is rather unlikely 
that we have missed any such long-period variable stars in For 4, if any is present. 

Both the Magellan and CTIO datasets were used, when available, to determine the 
pulsation characteristics of the variable stars from the study of their light curves. The 
possibility of combining the two datasets allowed us to increase the time-on-targets, greatly 
improving the phase coverage. In such a way we could solve the alias ambiguities at / ± 
1 cycle/day, which were very severe when considering only one dataset. The merged datasets 
also returned a better estimate of the period, epoch of maximum light and amplitudes of the 
light variation of the confirmed variable stars. We were able to derive reliable periods and 
light curves for 20 variable stars observed with both the Magellan and CTIO telescopes and 
for 3 stars with differential fluxes from the Magellan dataset only. For the remaining six cases 
we can only provide a rough estimate of the period or simply an identification of variability. 
The spectral window of our data is dominated by strong aliases at n cycles per year (with 
—3 < n < 3). We used period values to the fifth decimal place to fold the light curves, since 
this allowed us to significantly reduce the r.m.s. scatter of the truncated Fourier series best 
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fitting the data. However, our period determinations are not free from the annual aliasing 
effect, and this can easily affect the fourth decimal place. In other words, the true frequency 
could be an alias at n cycle/year of that corresponding to the periods finally adopted for the 
variable stars (see Table 3). The standard deviations of the least-square fits of the Magellan 
light curves are in the range from 0.02 to 0.04 mag and from 0.01 to 0.07 mag in V and B, 
respectively — which is quite acceptable and comparable with the expected precision of stars 
with V ~ 21.5 mag in a very crowded field. 

Twenty-seven of the variable stars detected in For 4 clearly appear to be RR Lyrae stars: 
22 ab-, 3 c-type, and 2 double-mode pulsators. However, among the listed fundamental-mode 
pulsators three have uncertain classifications, and a fourth one, star V13, has the typical 
period of a fundamental-mode pulsator, but the light curve shape of a first overtone variable 
star. Although the star shows significant gaps in the light curve, we could not fit the data 
with a shorter period. While long-period RRc stars have been observed in a few GC's 
(namely u Centauri = NGC 5139, NGC 6388, NGC 6441, and one star in M3; see Catelan 
2004b), such stars tend to be bright, whereas V13 appears to be even fainter than average, 
which is quite unusual and makes it very hard to explain its long period if it is indeed a 
first-overtone pulsator. Also, if [Fe/H] ~ —2, For 4 would by far be the lowest-metallicity 
globular cluster known to harbor a long-period RRc. Such circumstantial evidence may be 
in favor of our RRab, as opposed to a long-period RRc, classification for V13 (see Table 3). 
We note that an RRab star with behavior similar to V13 was observed in NGC 6441 by 
Pritzl et al. (2003). 

Coordinates for the variable stars are provided in Table 3, along with period, type, time 
of maximum light, number of phase points in the Magellan and CTIO datasets separately, 
intensity-averaged (V) and (B) magnitudes, and amplitudes of the light variation (A v and 
A B ). Examples of light curves are shown in Figure 2. The complete atlas of light curves is 
presented in the Appendix. Variable stars V8 and V9 have lights curve only in differential flux 
and large gaps around the maximum light (see Fig. 5). In turn, their period and classification 
in type are very uncertain. The folded light curves of V14 show scatter larger than for the 
other variable stars, in the CTIO dataset in particular (see Fig. 6), which renders the star's 
mean magnitudes and amplitudes uncertain. However, the individual measurements seem 
to confirm the intrinsic faintness of V14, which appears to be about 0.6 mag fainter than 
the average level of the cluster's HB. This faint magnitude would put V14 likely beyond the 
limits of Fornax if the star is a bona-fide RR Lyrae variable. Stars V10, V15 and V29 also 
are interesting given their small amplitudes for their periods. The Magellan light curves of 
V5 clearly show maxima at different heights, suggesting a double-mode nature, however we 
could not derive reliable periodicities for the star. 
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4.1. Membership Probability 

Given the different stellar content of For 4 and the surrounding Fornax field populations, 
estimating the expected number of RR Lyrae variable stars belonging to For 4 from the overall 
surface density profile is not straightforward (see Mackey & Gilmore 2003b). We have used 
both the outermost variable stars in the Magellan field and the variable stars found in the 
larger CTIO field to empirically estimate the expected surface density of RR Lyrae in the 
Fornax galaxy. In the Magellan data, 6 variable stars have been identified beyond 50", in 
an area of 1.23 x 10 4 arcsec 2 . On the basis of the overall surface density profile of For 4, we 
can assume that these variable stars are not associated with the cluster. Then we estimate a 
surface density of 4.9 x 10~ 4 RRL per arcsec 2 in the Fornax field. As for the CTIO data, RR 
Lyrae stars have been counted in an annulus between 140" and 220" from the center of For 4, 
which is an optimal compromise between the statistical benefits of sampling a large area and 
the need to keep our measurement local in view of the spatially varying stellar populations 
of Fornax. A surface density of 7.2 x 10~ 4 RRL arcsec -2 is inferred from the CTIO dataset. 
In the inner region, within a radius of 30" from the center, 1.4 and 2.0 field RR Lyrae stars 
are therefore expected, respectively. In conclusion, among the RR Lyrae stars within about 
30" from Fornax 4, about 2 stars possibly belong to the Fornax dSph field. 

The same exercise for the innermost variable stars found with ISIS shows that they are 
most likely to be all cluster members. 

5. Oosterhoff type and metallicity 

The pulsation properties of the For 4 RR Lyrae stars are summarized in Table 4, where 
we list the number of fundamental- mode (RRab), first overtone (RRc) and double-mode 
pulsators (RRd) in the first three columns, respectively; the ratio of number of RRc to total 
number of RR Lyrae stars (RR tot ), in column 4; the average periods of the ah- and c-type RR 
Lyrae stars, in columns 5 and 6; and the shortest and longest RRab periods, in columns 7 
and 8. The average period of the c-type RR Lyrae stars is (P c ) = 0.360 d (a — 0.042 d, 
average on 3 stars), or 0.374 d (a = 0.044 d, average on 4 stars), if the RRd star is included. 
The a6-type RR Lyrae stars have {P a b) — 0.594 d (a — 0.054 d, average on 19 stars). The 
quoted sigma values are the dispersion of the mean. V13 and V14 were not considered. 
Average and a do not change if V13 is included, but {Pab) = 0.596 d if V14 is added. The 
shortest-period RRab is star V16 with P a b,min = 0.5191 d, and the longest-period RRab is 
star V29 with P ab , max = 0.6986 d. 

The period distribution of the RR Lyrae stars in For 4 is shown in Figure 3. The 
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right panel shows basically the same distribution as in the left panel, only with the RRc 
periods "fundamentalized" according to the relation logPj = logP c + 0.128. The RRab 
distribution may give some hints of having two separate peaks at (P a b) = 0.55 d and (P a b) = 
0.64 d, respectively — thus resembling the case of NGC 1835 in the Large Magellanic Cloud 
(Soszynski et al. 2003). However, a KMM test (Ashman, Bird, & Zepf 1994) performed on 
our data shows that they are consistent with a unimodal distribution in P a b. Interestingly, 
the fundamentalized period distribution is peaked at a value 0.6 d, which is significantly 
longer than for either Ool or OoII clusters. We refer the reader to Catelan (2004a) for an 
extensive discussion of the implication of the peaked distributions in fundamentalized RR 
Lyrae periods. 

According to the RR Lyrae star average periods For 4 appears to be an Oosterhoff- 
intermediate (Oo-Int) cluster, an Oo-type that has very few counterparts among the Milky 
Way GC's (Catelan 2004b, 2005). Even if we only consider stars in the inner region of the 
cluster (R < 30") and if we also take into account that one or two of them could not belong 
to the cluster (see Section 4.1), we have (P a b) = 0.594 d, thus confirming a true Oo-Int 
behavior, similarly to what was found for the other 4 GC's in Fornax by Mackey & Gilmore 
(2003a). 

Period-amplitude distributions for the For 4 RR Lyrae stars with light curves in magni- 
tude units are shown in Figure 4. Expanded symbols are used for stars located at distances 
less than 30" from the cluster center. Figure 4 shows that the majority of the a6-type RR 
Lyrae stars in For 4 appears to be on the Ool line. It is also noteworthy that the ratio of 
RRc's to the total number of RR Lyrae stars is more like Ool (see Table 4). In Figure 4 
there are, however, a number of RRab's that pile up near the OoII location. However, the 
two stars closer to the OoII line are located outside R = 30", thus their membership is less 
certain. 

As discussed in the Introduction, the metal abundance of For 4 is somewhat uncertain, 
being [Fe/H] = —2.0 in the Buonanno et al. (1999) study, and [Fe/H] = —1.5 in Strader 
et al. (2003). There are some metallicity estimates that can be made using the pulsation 
properties (namely periods and amplitudes) and the shape of the light curve (namely the 
parameters of the Fourier decomposition of the light curve) of the RR Lyrae stars. Applying 
the Alcock et al. (2000) relation, which is based on the period and V-band amplitude of 
RRab stars, we found a mean metallicity of —1.68 for the cluster (or —1.66 for the 10 RRab 
stars within 30" from the cluster center). This is intermediate between the Buonanno et 
al. and Strader et al. values. Using relations from Sandage (2006), we find [Fe/H] values 
of —1.96 from the mean period of the fundamental-mode RR Lyrae stars, —2.25 from the 
shortest-period RRab (V16, a star located within 30" from the cluster center, hence very 
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likely belonging to For 4), and —1.38 from the longest-period RRab (star V29, the most 
external of the variable stars in our sample, hence probably belonging to the Fornax dSph 
field). Finally, from Sandage (1993), we find [Fe/H] = —1.90 from the mean period of the 
three cluster RRc stars, which are all located within 30" from the For 4 center. 

According to Jurcsik & Kovacs (1996) the [Fe/H] abundance of a fundamental-mode 
RR Lyrae star is a linear function of the period P and of the parameter 31 = 3 — 30i of 
the Fourier decomposition of the star's l^-band light curve. Individual metal abundances 
were derived for the For 4 fundamental mode RR Lyrae stars by applying the Jurcsik & 
Kovacs method (equation 3 in their paper) to the RRab variable stars whose light curves 
satisfy the compatibility conditions (namely light curve completeness and regularity criteria) 
for the Fourier parameters to predict reliable empirical quantities (see Kovacs & Kanbur 
1998)@ Two of the For 4 RRab stars are found to satisfy these conditions, namely star V18 
and V24, for which we derive [Fe/H] = —1.62 and —1.82, respectively. These photometric 
metallicities are on the Jurcsik (1995) metallicity scale and can be transformed to the Zinn & 
West (1984) scale by applying the transformation relation provided by equation 4 in Jurcsik 
(1995), leading to [Fe/H] = —1.75 and —1.89, respectively. 

Thus, different methods based on the cluster RR Lyrae stars seem to suggest a low 
metal abundance for For 4, closer to the Buonanno et al. (1999) estimate. However, we 
should be aware that the calibrations of both the period-amplitude- [Fe/H] (P — A — [Fe/H]) 
and the period-03i-[Fe/H] (P — </> 31 — [Fe/H]) relations are based on MW GC's that do follow 
the Oosterhoff dichotomy, and might thus not be suitable to derive the metal abundance 
of an object that does not follow the dichotomy such as For 4. Also, the P — A — [Fe/H] 
relations have a known scatter at any given [Fe/H]. For instance, there is a range of [Fe/H] 
where Ool and OoII objects overlap (at around the M2/M3/M5/M62 metallicity), and there 
is accordingly a tendency to assign a much lower metallicity for an OoII or Oo-Int object in 
this domain than should really be the case. 

Nevertheless, we note that none of the conclusions on the Oo-Int status of For 4 are 
affected in any way by a change in the cluster metal abundance from —2.0 to —1.5 dex. 
Indeed, for both values of metallicity For 4 is found to fall in the Oo-Int region in the 
(Pab) — [Fe/H] diagram (Pritzl et al. 2002; Catelan 2005). In particular, adopting a metallicity 



2 According to Kovacs & Kanbur (1998), deviations of the Fourier parameters should not exceed the 
maximum value (D m ) of 3, with maximum deviations D m > 3 possibly indicating incompatibility with the 
empirical predictions. The deviation parameters Dp are defined as Dp —\ F t, s — F ca i c | /o\f, where F i, s , 
Fcaic are respectively the observed value of a given Fourier parameter and its predicted value from the other 
observed parameters, and ap is the respective standard deviation (see eq. 6 and Table 6 of Jurcsik & Kovacs 
1996). 
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of [Fe/H] = — 1.5 according to Strader et al. (2003), would place For 4 right in the middle 
of the Oo-Int band in the [Fe/H]-HB type diagram (see Fig. 8 in Catelan 2005), whereas 
adopting [Fe/H] = —2.0 actually places the cluster slightly to the left of this band. In 
fact, and as discussed in a companion paper (Catelan et al. 2007, in preparation), the 
key quantities defining Oosterhoff status appear to be the average and minimum periods 
of the a6-type RR Lyrae stars. In terms of both these quantities, For 4 appears entirely 
consistent with an Oo-Int classification. This is also what is generally expected on the 
basis of the position of the cluster in the HB type-metallicity plane. Note, in this sense, 
that our photometry and preliminary number counts suggest a revised Lee-Zinn HB type of 
(B — R)/(B+V+ R) —0.75 for the cluster — where B, R, and V represent the numbers of 
blue, red, and variable stars on the HB, respectively. The Oo-Int behavior further confirms 
the broad similarity of For 4 to the MW cluster Rup 106. Indeed, assuming (P a b) — 0.616 d 
(Kaluzny, Krzeminski, & Mazur 1995) and [Fe/H] = —1.67 (Harris 1996), Rup 106 is one 
of the few Galactic GC's falling in the Oosterhoff-intermediate band in the [Fe/H]-HB type 
diagram. However, there are some differences between the two clusters at a more detailed 
level: for instance, Rup 106 has no RRc whatsoever, while For 4 has a few; and, if we adopt 
a boundary of 0.62 d between Oo-Int and OoII clusters (Catelan 2004b, 2005), Rup 106 is 
in fact a "borderline" Oo-Int cluster, whereas For 4 falls closer to the middle of the Oo-Int 
zone. 



6. Cluster distance 



The average apparent magnitude of the RR Lyrae stars within 30" from the For 4 
center (and hence more likely cluster members), and excluding star V14 which may be a 
background object, is (V(RR)) = 21.43 ± 0.03 mag (a = 0.10 mag, average over 12 stars) 
and (B(RR)) = 21.76 ± 0.04 mag (er = 0.13 mag, average over 12 stars). Average values 
become instead (V(RR)) = 21.38 ± 0.02 mag (a = 0.10 mag, average over 19 stars) and 
(B(RR)) = 21.74 ± 0.03 mag (er = 0.11 mag, average over 19 stars) if we consider "all" the 
variable stars in the field of For 4 but V14. Both these values, and the latter in particular, 
are consistent with the value of (V(HB)) = 21.37±0.04 mag, for the mean level of the (red) 
old HB stars in the Fornax field, as derived by Saviane, Held, & Bertelli (2000). On the other 
hand, our values are respectively 0.09 and 0.14 mag brighter than Vhb = 21.52 ± 0.05 mag 
found by Buonanno et al. (1999) using HB stars within a distance of 18" from the For 4 
center. The 0.09 mag difference remains even if we consider only the 8 RR Lyrae stars 
with reliable photometry, located within 18" from the cluster center. Unfortunately, we lack 
reliable photometry for a further 9, more internal, RR Lyrae stars which are within 8" from 
the cluster center. We note that a ~ 0.1 mag difference between our (V(RR)) and Buonanno 
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et al.'s Vhb might easily be accounted for if the For 4 RR Lyrae stars were more evolved 
than their non-variable HB counterparts, as observed in a number of the metal-poor GC's 
in the MW. Unfortunately, unlike those GC's, For 4 has a predominantly red HB, although 
it should be noted that the Catelan (1993) synthetic HB's for a predominantly red HB 
population also predict a difference in luminosity, at a fixed temperature, between the lower 
envelope of the RR Lyrae distribution on a CMD and the zero-age HB luminosity level (see 
his Table 1, last row, and also Fig. 3 in Catelan 1992). 

In order to derive the distance to For 4 from the average luminosity of its RR Lyrae 
stars we need an estimate of the cluster reddening and metallicity, as well as values for the 
slope and zero point of the RR Lyrae luminosity-metallicity relation. We have assumed an 
absolute magnitude of My = 0.59 ± 0.03 for RR Lyrae stars of [Fe/H] = —1.5 (Cacciari 
& Clementini 2003; Catelan 2005), and AM y /A[Fe/H] = 0.214 ± 0.047 (Clementini et al. 
2003, Gratton et al. 2004) for the slope of the luminosity-metallicity relation. 

We then adopt a standard extinction law [A v = 3.1 x E(B — V)] and a reddening value 
of E(B-V) = 0.10 ± 0.02 for the cluster, which is the weighted average of E(B-V) = 
0.12 ± 0.05 from Mackey & Gilmore (2003a) and E(B-V) = 0.08 ± 0.03, the reddening we 
find by matching the blue edge of the RR Lyrae strip in M3 [E(B — V) M 3 = 0.01; Harris 
1996] to the colors of the bluest RR Lyrae stars in For 4. 

For [Fe/H] = —2.0 (Buonanno et al. 1999), the cluster's distance modulus is then 
20.64 ±0.09 mag, while it becomes 20.53 ±0.09 mag if the Strader et al. (2003) metallicity is 
adopted. Here, errors are the sum in quadrature of uncertainties of 0.03 mag in (V(RR)) (dis- 
persion of the average), 0.05 mag in the zero point of the photometry, 0.02 mag in E(B — V) 
(corresponding to 0.06 mag in Ay), and of 0.03 mag and 0.047 mag/dex, respectively, in the 
zero point and in the slope of the RR Lyrae luminosity-metallicity relation. 

Previous distance estimates for the Fornax dSph field range from /x = 20.59 ±0.22 mag 
(Buonanno et al. 1985) to 20.76 mag (Demers et al. 1990; Buonanno et al. 1999, for an 
assumed metal abundance of the Fornax field stars of —1.4). Saviane et al (2000) infer fiQ = 
20.70 ± 0.12 mag from the tip of the field stars' red giant branch and fx = 20.76 ± 0.04 mag 
from the mean magnitude of Fornax's old HB stars, which are assumed to have [Fe/H] ~ —1.8 
and E(B — V) = 0.03 mag. Buonanno et al. (1999) derive an average distance modulus of 
/i = 20.62 ± 0.08 mag for the Fornax clusters #1, 2, 3, and 5 and, more recently, Mackey & 
Gilmore (2003a) find for the same clusters distance moduli in the range from 20.58 ±0.05 to 
20.74 ± 0.05 mag, with an average value of /x = 20.66 ± 0.03 (random) ±0.15 (systematic) 
mag. 

For ease of comparison we have summarize in Table 5 the various distance determina- 
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tions and the different assumptions about reddening, metallicity and RR Lyrae luminosity- 
metallicity relation they are based on. Our "long" distance modulus of 20.64 ± 0.09 mag, 
based on [Fe/H] = —2.0 and E(B — V) = 0.10 mag for For 4, agrees well with both the 
Saviane et al. (2000) and the Mackey & Gilmore (2003a) estimates, once differences in 
the adopted reddening, cluster metallicity, and RR Lyrae luminosity-metallicity relation are 
properly accounted for. It is, however, 0.09 mag shorter than one would infer by assuming 
the Buonanno et al. (1999) V HB = 21.52 ± 0.05 mag for the cluster. On the other hand, 
the "short" distance modulus of 20.53 ± 0.09 mag, based on [Fe/H] = —1.5 for For 4, is 
significantly shorter than any previous distance determination for the Fornax system. 



7. Summary and Conclusions 

We have presented the first study of the variable star population in For 4, the GC 
located in the central region of the Fornax dSph galaxy. Our sample includes 29 variable 
stars, of which 27 are RR Lyrae stars (22 fundamental-mode, 3 first overtone, and 2 double- 
mode pulsators). Twenty- two of these variable stars are located within 30" from the cluster 
center, hence are very likely cluster members. The average periods of ab- and c-type RR 
Lyrae stars, and the minimum period of the a6-type pulsators, point to an Oo-Int status for 
the cluster, unlike what is seen for the vast majority of the Galactic GC's. A similar Oo-Int 
classification is suggested by Mackey & Gilmore (2003a) also for the other 4 GC's in Fornax. 

The Oo-Int type of For 4 further supports the cluster similarity with the unusual "young 
halo" Galactic cluster Rup 106, and strengthens the claim for an extragalactic origin of this 
MW cluster. Our results on the Fornax field and cluster variable stars also indicate that 
while galaxies like the Fornax dSph unlikely contributed large fractions of the Galactic halo, 
they may have provided instead some of the unusual clusters being observed in the MW (see 
e.g. Catelan 2007). 

The average luminosity of the RR Lyrae stars that are more likely cluster members 
is (V(RR)) = 21.43 ± 0.03 mag, leading to distance moduli of /i = 20.64 ± 0.09 mag or 
fiQ = 20.53 ± 0.09 mag, depending on whether Buonanno et al.'s (1999) [Fe/H] = —2.0 or 
Strader et al.'s (2003) [Fe/H] = —1.5 determination is adopted for the cluster metallicity, 
and for an assumed reddening E(B — V) = 0.10 mag. The properties of the cluster RR Lyrae 
stars also support a low metal abundance for For 4 closer to Buonanno et al. value, and 
a distance modulus of /io = 20.64 ± 0.09 mag (D = 134 ± 6 kpc). This distance modulus 
agrees well with both Saviane et al. (2000) and Mackey & Gilmore (2003a), once differences 
in the adopted reddening, cluster metallicity, and adopted RR Lyrae luminosity-metallicity 
relation are properly taken into account and corrected for. 
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Mackey & Gilmore (2003a) find that the distance moduli of the Fornax clusters #1, 
2, 3 and 5 are consistent with a line of sight depth of ~ 8 — 10 kpc for this galaxy. The 
distance modulus of For 4 and our preliminary results from the study of the RR Lyrae stars 
in clusters #2, 3 and 5 show that For 4 is at the center of the Fornax GC's distance moduli 
distribution, and confirm Mackey & Gilmore's finding of a line-of-sight depth in Fornax, 
although reducing its extent to about 7 kpc. Further insight on the line of sight depth of 
the Fornax dSph will be gained from our study of the galaxy's field variable stars. 
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A. Atlas of the Light Curves 

Atlas of light curves for the 29 variable stars identified in For 4. Photometric data 
are folded according to the ephemerides provided in Table 3. Variable stars are ordered 
by increasing distance from the cluster center. Only differential flux light curves from the 
Magellan V dataset are available for stars from VI to V9. 
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Table 1: Instrumental set-ups and logs of the observations 



Dates 


Telescope 


Instrument 


Detector 




Resolution 




FOV 


N B 


N v 


Photometric precision 






















(HB level) 


UT 






(pixel) 




(///pixel) 










(mag) 


Nov. 13 - 15, 2003 


Magellan/Clay 


MagIC 


2048 X 2048 


SITc 


0.069 


2.35' 


' X 2.35' 


10 


39 


0.01 


Dcc.l - 2, 2004 


Magellan/ Clay 


MagIC 


2048 X 2048 


SITc 


0.069 


2.35' 


' X 2.35' 


9 


19 


0.01 


Nov. 14, 2003 


Blanco/CTIO 


MosaicII 


2048 X 4096 


SITc 


0.27 


36' 


X 36' 


8 


16 


0.03 - 0.05 


Oct. 25 - 26, 2004 


Blanco/CTIO 


MosaicII 


2048 X 4096 


SITc 


0.27 


36' 


X 36' 


51 


124 


0.03 - 0.05 


Sept. 9, 2005 


Blanco/CTIO 


MosaicII 


2048 X 4096 


SITc 


0.27 


36' 


X 36' 


2 


2 


0.03 - 0.05 



Table 2: V, B photometry of the For 4 variable stars with light curves in magnitude scale 



Star V18 - RRab 



HJD 


V 


HJD 


B 


(-2452956) 






(-2452956) 




0.697471 


21 


.86 


0.749611 


22.26 


0.702402 


21 


.79 


1.563541 


21.66 


0.708871 


21 


.76 


1.604275 


21.93 


0.714947 


21 


.79 


1.631649 


21.97 


0.729866 


21 


.77 


1.665676 


22.10 


0.735942 


21 


.83 


1.690539 


22.18 


0.742018 


21 


.83 


1.717479 


22.24 


0.757122 


21 


.84 


1.744663 


22.33 


0.763198 


21 


.84 


1.772030 


22.31 


0.784262 


21 


.84 


2.554781 


20.88 



A portion of Table 2 is shown here for guidance regarding its form and content. The entire 
catalog is available in the electronic edition of the journal. 
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Table 3: Identification and properties of the For 4 variable stars 



Name 
(a) 


Id 




(2000) 




S 

(2000) 




Type 


P 

(days) 


Epoch 
(-2450000) 


<v> 


N v 


<B> 


N B 




Av 


A B 


Notes 


VI 


LC378 


2 


40 : 


07 


62 


-34 


: 32 : 


09 


4 


RRab 





54 :: 








54 


















(b,c) 


V2 


LC360 


2 


40 : 


07 


70 


-34 


: 32 : 


09 


8 


RRab 





6261 


2955.77 






52 


















(b) 


V3 


LC330 


2 


40 : 


07 


47 


-34 


: 32 : 


10 


7 


RRab 





5968 


3342.640 






56 


















(b,c) 


V4 


LC386 


2 


40 : 


07 


75 


-34 


: 32 : 


09 





RRab 





65 :: 








53 


















(b,c) 


V5 


LC319 


2 


40 : 


07 


80 


-34 


: 32 : 


11 





RRd? 





46 :: 








57 


















(b) 


V6 


LC381 


2 


40 : 


07 


33 


-34 


: 32 : 


09 


2 


RRab 





52936 


2956.736 






55 


















(b) 


V7 


LC178 


2 


40 : 


07 


49 


-34 


: 32 : 


15 


1 


RRab 





5 : : 








55 


















(b,c) 


V8 


LC486 


2 


40 : 


07 


40 


-34 


: 32 : 


05 


1 


RR? 





83 :: 








57 


















(b,c) 


V9 


LC291 


2 


40 : 


07 


99 


-34 


: 32 : 


11 


5 


RR? 





69 :: 








57 


















(b,c) 


V10 


2658 


2 


40 : 


07 


34 


-34 


: 32 : 


03 


6 


RRab 





64450 


2956.740 


21 


31 


51 




21 


75 


19 







40 


0.49 




Vll 


2051 


2 


40 : 


08 


39 


-34 


: 32 : 


09 


7 


RRab 





53085 


3342.568 


21 


42 


57 


+ 136 


21 


80 


17 + 


60 


1 


06 


1.36 




V12 


2272 


2 


40 : 


07 


99 


-34 


: 32 : 


20 


5 


RRab 





59506 


3342.545 


21 


35 


58 


+ 139 


21 


74 


19 + 


61 





83 


1.03 




V13 


2949 


2 


40 : 


06 


81 


-34 


: 32 : 


16 


9 


RRab 





58767 


2958.620 


21 


57 


52 




21 


79 


13 







85 


0.75 :: 


(d) 


V14 


1956 


2 


40 : 


08 


55 


-34 


: 32 : 


12 


4 


RRab 





64520 


3303.680 


22 


04 


55 


+ 109 


22 


27 


19 + 


47 


1 


14 


1.40 




V15 


2076 


2 


40 : 


08 


33 


-34 


: 32 : 


19 


8 


RRab 





67000 


2958.640 


21 


45 


58 


+ 139 


21 


96 


18 + 


58 





27 


0.42 




V16 


2651 


2 


40 : 


07 


35 


-34 


: 31 : 


55 


7 


RRab 





51910 


3341.555 


21 


49 


57 




21 


80 


17 







64 


0.83 




V17 


2664 


2 


40 : 


07 


33 


-34 


: 32 : 


25 


9 


RRc 





31618 


2958.573 


21 


30 


57 


+ 141 


21 


55 


17 + 


63 





53 


0.67 




V18 


3177 


2 


40 : 


06 


38 


-34 


: 32 : 


15 


7 


RRab 





56297 


2958.5475 


21 


44 


58 


+ 124 


21 


79 


19 + 


39 


1 


15 


1.54 




V19 


3303 


2 


40 : 


06 


14 


-34 


: 32 : 


08 


8 


RRab 





52651 


2958.600 


21 


55 


58 


+ 136 


21 


90 


19 + 


52 


1 


20 


1.50 




V20 


1685 


2 


40 : 


09 


09 


-34 


: 32 : 


17 


2 


RRc 





36479 


3342.544 


21 


54 


52 


+ 141 


21 


87 


18 + 


51 





40 


0.43 




V21 


2948 


2 


40 : 


06 


80 


-34 


: 32 : 


33 


9 


RRc 





39962 


3341.663 


21 


30 


56 


+ 131 


21 


62 


19 + 


50 





36 


0.42 




V22 


3544 


2 


40 : 


05 


68 


-34 


: 32 : 


20 


4 


RRab 





65255 


3341.540 


21 


42 


57 




21 


83 


19 + 


49 





57 


0.71 




V23 


1130 


2 


40 : 


10 


29 


-34 


: 32 : 


12 


4 


RRab 





55746 


3342.600 


21 


28 


56 


+ 136 


21 


76 


19 + 


56 





54 


0.69 




V24 


4153 


2 


40 : 


04 


39 


-34 


: 32 : 


40 


9 


RRab 





61631 


2958.695 


21 


30 


58 


+ 140 


21 


66 


19 + 


59 





69 


0.90 




V25 


2219 


2 


40 : 


07 


91 


-34 


: 33 : 


07 


2 


RRd 





41522 


2958.570 


21 


31 


56 


+ 137 


21 


64 


19 + 


61 





31 


0.56 


(e) 


V26 


4779 


2 


40 : 


03 


13 


-34 


: 31 : 


43 


4 


RRab 





61610 


2958.600 


21 


31 


57 




21 


66 


19 




1 


18 


1.46 




V27 


4989 


2 


40 : 


02 


63 


-34 


: 32 : 


42 





RRab 





59643 


3341.620 


21 


36 


58 


+ 140 


21 


60 


19 + 


61 





84 


1.12 




V28 


437 


2 


40 : 


11 


81 


-34 


: 31 : 


13 


6 


RRab 





55954 


3303.880 


21 


26 


55 


+ 119 


21 


62 


17 + 


51 





97 


1.24 




V29 


5361 


2 


40 : 


01 


81 


-34 


: 32 : 


54 


7 


RRab 





69860 


2958.587 


21 


23 


58 


+ 134 


21 


64 


19 + 


61 





40 


0.65 





a Variable stars were assigned increasing numbers starting from the cluster center that was set at 
a = 02:40:07.6, 5 = -34:32:10.0 (J2000). Stars from VI to V22 are located within 30" from the 
For 4 center, hence are more likely cluster members. 
b Light curves available in differential flux only. 
c Stars with gaps in the folded light curve. 

d Star with the typical period of a fundamental-mode pulsator, but with the light curve of a first 
overtone pulsator. 

e Double-mode RR Lyrae star with fundamental-mode period Pq = 0.56284 d and period ratio 
Pi/Pq = 0.7377; in the table we list the star's first overtone period. 
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Table 4: Average quantities for the For 4 RR Lyrae stars 



NRRab 


N RRc 


N RRd 


N RRc /N(RR tot ) 


(Pab) (d) 


(P c ) (d) P ab , min (d) P ab , max (d) 


(a) 


(a) 


(a) 


(a,b) 


(c) 


(d) 


22(16) 


3(3) 


2(1) 


0.12(0.16) 


0.594 


0.360 0.5191 0.6986 








0.19(0.20) 


0.596 


0.374 



a Values in parentheses consider only RR Lyrae stars within 30" from the cluster center. 

b Ratio of RRc to total number of RR Lyrae stars with (2nd raw values) and without (1st raw 

values) the RRd stars. 

c Fundamental-mode average period with (2nd raw value) and without (1st raw value) stars V13 
and V14. 

d First overtone average period with (2nd raw value) and without (1st raw value) the double mode 
star V25. 



Table 5: Distance determinations based on the HB luminosity of the Fornax dSph field and 



globular clusters 


V HB 


E(B-V) 


[Fe/H] 


M v vs. [Fc/H] 


M0 


Reference 


(mag) 


(mag) 










(field) 


0.03 






20.59±0.22 


Buonanno ct al.(1985) 


21.37±0.04 (field) 


0.03 


-1.8 


0.17x([Fc/H]+1.5)±0.57 


20.76±0.04 


Savianc ct al. (2000) 


(field) 


0.05 






20.76 


Demers et al.(1990) 


(field) 


0.05 




0.17x([Fc/H]+1.5)±0.57 


20.76±0.10 


Buonanno et al.(1999) 


21.28±0.01 (field) 


0.04±0.03 


-1.77 ±0.20 


0.22x([Fe/H]+1.5)+0.50 


20.72±0.10 


Greco et al. (2007) 


21.20±0.05 (For 1) 


0.04±0.05 


-2.20 ±0.20 


0.17x([Fc/H]+1.5)±0.57 


20.62±0.08 


Buonanno et al.(1999) 


21.35±0.05 (For 2) 


0.07±0.05 


-1.78 ±0.20 






Buonanno et al.(1999) 


21.20±0.05 (For 3) 


0.04±0.05 


-1.96 ±0.20 






Buonanno et al.(1999) 


21.20±0.05 (For 5) 


0.06±0.05 


-2.20 ±0.20 






Buonanno ct al.(1999) 


21.27±0.01 (For 1) 


0.07±0.01 


-2.05 ±0.10 


0.23x([Fc/H]+1.6)±0.56 


20.58±0.05 


Mackey & Gilmorc (2003a) 


21.34±0.01 (For 2) 


0.05±0.01 


-1.83 ±0.07 




20.67±0.05 


Mackey & Gilmore (2003a) 


21.24±0.01 (For 3) 


0.04±0.01 


-2.04 ±0.07 




20.66±0.05 


Mackey & Gilmore (2003a) 


21.33±0.01 (For 5) 


0.03±0.01 


-1.90 ±0.06 




20.74±0.05 


Mackey & Gilmorc (2003a) 


21.21±0.02 (For 3) 


0.04±0.03 


-1.91 ±0.20 


0.22x([Fe/H]+1.5)+0.50 


20.68±0.11 


Greco et al. (2007) 


21.43±0.03 (For 4) 


0.10±0.02 


-2.01 ±0.20 


0.214x([Fc/H]+1.5)±0.59 


20.64±0.09 


This paper 
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Fig. 1. — Color- magnitude diagram of For 4 from the Magellan data. Variable stars are 
plotted according to their intensity-averaged magnitudes and colors, using different symbols 
for the different types. Open circles: ab-type RR Lyrae stars; filled circles: first overtone RR 
Lyrae stars; asterisk: double-mode pulsator (RRd). The open triangle marks a peculiar RR 
Lyrae star (V13). 
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Fig. 2. — V, B light curves of RR Lyrae stars in For 4. Left panels, from top to bottom: an 
RRc, an RRab and an RRd star from the Magellan dataset only. Right panels: the same, 
for the combined Magellan + CTIO dataset. 
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Fig. 3. — Left panel: Period distribution of fundamental- mode and first overtone RR Lyrae 
stars in For 4. Right panel: Fundamentalized period distribution for RR Lyrae stars in For 4. 
Interestingly, the peak of the distribution, at 0.6 d, is seemingly located at a longer period 
than for either Ool or OoII GC's. 
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log P 



Fig. 4. — V, B period-amplitude diagrams of For 4 RR Lyrae stars. Symbols are as in 
Figure 1, but those indicating variable stars located at distances less than 30" from the 
cluster center are expanded. The straight lines are the positions of the Ool and OoII Galactic 
GC's according to Clement & Rowe (2000). Period-amplitude distributions of the bona 
fide regular (solid curves) and well-evolved (dashed curves) ah RR Lyrae stars in M3 from 
Cacciari, Corwin, & Carney (2005) are also shown for comparison. 



-28- 




Fig. 5. — Light curves in differential flux, from the Magellan V dataset. 
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Fig. 6. — V (left panels), B (right panels) light curves from the combined Magellan {solid 
circles) plus CTIO {open circles) datasets. 
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Fig. 6. — continued 
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Fig. 6. — continued 
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